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ABSTRACT 

We report results of an observation of the 7-s pulsar 4U1 626-67 
with the A2 experiment on HEAO 1. The phase-averaged X-ray spectrum 
between 0.7 and 60 keV is complex as are the constituent spectra which 
change radically as a function of pulse phase. Included in this spectral 
change is the sudden appearance (A<f> <.0.1) and subsequent decay of a 
continuum or emission feature with a mean energy of 19 keV which 
contains about h the power in this spectral range. Pulse timing results 
include a new determination of the pulse period and a factor 8 reduction 
in the upper limit for the light travel time for orbital periods between 
1 and 7 hours. We discuss these findings for this system and compare 
them with the general nature of pulsar spectra. 
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I. INTRODUCTION 

It has been established that X-ray pulsars belong to a spectral 
class which is characterized by an extremely hard spectrum (power law 
index a £ 1) in the 2-15 keV energy band (e.g. Ulmer 1975). The 7-sec 
pulse period of 4U1626-67 was discovered when a crude study of its spectrum 
(Markert et al. 1977) indicated that it was a member of this spectral 
class and a good candidate for a period search (Rappaport et al. 1977). 

If this pulsar Is a member of a binary system, then variations in the 
pulse arrival times due to radial orbital motion can determine the orbital 
light travel time (a x sini), the orbital period (P Qrb )» and hence the 
mass function for the system. Rappaport et al. did not detect this effect 
but were able to exclude a significant amount of the a x sini vs. P Qrb 
plane in parameter space. This result and that of McClintock et al. 

(1977) which suggested a faint blue star as the optical counterpart, led 
to a model of a binary system containing a neutron star and a low mass dwarf 
companion with P Qrb £ 7 hours (Joss, Avni, and Rappaport 1978). The 
optical identification was confirmed by Ilovaisky, Motch, and Chevalier 
(1978) when they detected 7-sec pulsations from the star. 

The pulse profile is very energy dependent showing drastic changes 
in amplitude, phase, and structure over a relatively small energy range 
(Rappaport et al. 1977). However there have been to date no detailed 
spectral studies of this source. We present here a HEAO 1 pointed 
observation of 4U1 626-67 from which we have been able to obtain pulse 
phase spectroscopy and further improve the upper limits to orbital motion. 


— 3 - 

II. THE OBSERVATIONS 

The HEAO 1 A2 + X-ray spectroscopy experiment (Rothschild et al. 1978) 

+ The A2 experiment oh HEAO-1 is a collaborative effort led by E. Boldt of 
GSFC and 6. Garmire of CIT, with collaborators at GSFC, CIT, JPL and UCB. 

observed 4U1 626-67 for 5 hours starting at 15 hrs UT on 1978, March 
29. Three co-aligned proportional counters were utilized, each containing 
a different filling gas. They were 1) propane filled (LED 0.7 to 1.9 keV at 
this time), 2) argon filled (MED 1.7-18 keV), and 3) xenon filled (HED 3.- 
80. keV). The MED and HED counts were pulse height analyzed into 64 
channels with a temporal resolution of 20 ms. The broad band accumulation 
interval of the LED was 80 ms with 10 channel pulse height analysis every 
40 s. The diffuse and particle backgrounds were obtained from off source 
observations prior to this measurement. The LED background was determined 
using a field-of-view subtraction technique. 

III. THE PULSE PERIOD 

Periodogram analysis of the data revealed the 7.6 s pulse to be clearly 
present in each of the detectors. The energy dependent nature of the 
pulse profile was also evident. In order to maximize the amount of pulse 
phase information available the counts were binned into three energy ranges: 
0.7 to 1.9 keV; 3. to 14. keV; and 14 to 30 keV. Each 30 minutes of data 
were folded about the period determined by the periodogram analysis. 

The phase of each cf these light curves relative to an arbitrary reference 
point was then determined by fitting each light curve to a standard template. 
This process was repeated twice, the first time with one of the 30 min 
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folds used as the template, the second time using a master template con- 
structed by superposing each of the 30 min light curves (excluding the 
light curve being fitted) shifted according to the phases determined from 
the first Iteration. The light curves constructed from the second 
Iteration are shown in Figure 1. Six measurements of the pulse phase 
were obtained In each energy Interval with one sigma uncertainties of 
0.07 It s In both the 3-14 and 14-30 keV bands and 0.15 It s In the 0.7 
to 1.9 keV band. 

Joss et al. (1978) using SAS 3 set stringent upper limits to orbital 
motion with a period greater than 7 hrs. However for periods less than 
this their upper limit to the projected semi-major axis (a x slnl) was 
of the order of only 0.8 It s (Figure 2). The gross effective area of 
the A-2 detectors is * 20 times greater than SAS. By combining the phase 
measurements from each energy band we see no evidence for orbital periods 
less than 7 hrs with a x sini K 0.1 It s/(2 sigma) (Fig. 2). 

The heliocentric pulse period derived from these measurements was 
7.678174 + 0.000051 (2 sigma). Joss et al. obtained a value of 7.6806273 
+ 0.0000005 a year earlier, which yields a P/P of -3.2 x 10” 4 yr" 1 . This 
Is almost twice that measured by Joss et al. but variations of this 
magnitude are typical in the spinup rates of other binary pulsars 
(e.g. Mason 1977). 

IV. AVERAGE SOURCE PROPERTIES 

The average X-ray luminosity of 4U1626-67 during the pointed obser- 
vation was 2.9 (+ 0.3) x 10 33 d^ pc erg sec"^ between 0.7 - 60 keV, This 
result Is similar to that observed with UHURU (Giacconi et al . 1974) and 
OSQ-7 (Markert et al. 1977) in the 2-6 keV range, and also consistent with 
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the result of Rappaport et al. (1977) with the 1-*C keV energy range of 
SAS 3. For one week In September 1 977 we obtained scanning observations of 
4U1626-67. The average intensity was a third lower to that observed six 
months later and a factor * 2 variability was observed on timescales of 
hours. We can not comment on the existence of the * 10 sec quasi-periodic 
intensity oscillations discussed by Joss et al. because systematic aspect 
uncertainties during HEAO 1 pointings can result in spurious effects 
with a similar time scale. The scanning data consists of * 2 minutes of 

3 

observation every * 35 minutes and is also not relevant for the * 10 
sec timescale. 

The phase-averaged spectrum between 0.7-60 keV is shown in Figure 3. 
Fravdo et al. (1978) discuss the method of spectral analysis. This spectrum 
is quite complex and can be fitted by a range of different models. 

Between 5-20 keV a power luw number index of * 0.4 adequately describes 
the data. Above 20 keV the spectrum falls off sharply. This is a canonical 
phase-averaged pulsar spectrum. Below 5 keV an additional thermal or 
blackbody continuum component is present in some models. The apparent 
column density of cool material in the line of sight (Brown and Gould 
1970; Fireman 1974) is model -dependent but is always £ 4 x 10 H atoms 
cm“ 2 . A very broad line near 6.5 keV (FWHM K 2 keV) with .001 -.01 photons 
cm sec significantly improves the model fits (ax£ * 1) and could 
indicate the presence of broadened iron line emission which has been 
seen in other pulsars (see Pravdo 1979 and references therein). However, 
the average spectrum is composed of radically varying constituent spectra, 
and thus the models for it may be misleading. 
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V. PULSE PHASE SPECTRA OF 4U1 626-67 

Although the broad-band Intensity pulses shown In Figure 1 are a 
representation of the spectral evolution, the pulse features (e.g. the rise 
time) are sensitively dependent on the chosen broad-band energy limits. 

In a mode! In which the X-ray emission Is angle-dependent (the super- 
strong magnetic field model e.g. Basko et al. 1975) there is more 
Information In the energy distribution of photons at a particular phase 
than with the distribution of photons at a particular energy over all 
phases. Interpretation of the latter Is difficult since one can not dis- 
tinguish the separate effects of energy- Independent intensity changes 
(an overall pulse envelope) and spectral changes with pulse phase. We 
show in Figure 4 the evolution of the 4U1626-67 spectrum across the pulse. 
Each spectrum, progressing down the left column and then the right column. 

Is comprised of 0.1 phase. There is bimodal distribution among the 10 
spectra displayed. Figure 5 shows spectra 3 and 8 superposed. Spectra 1-5 
and 10 harden considerably with the addition of a continuum or broad emission 
feature centered near 19 keV. There is also a relative deficiency near 10 
keV in these spectra. Again, as in the case of the phase-averaged spectrum, 
no simple continuum model can fit these data. We present examples of two 
possible classes of models. 

A superposition of continuum components plus a high energy cutoff 
can describe the spectra in Figure 4. Table 1 lists relevant parameters 
(a-power law index, E c and Ep-cutoff parameters, T-temperature and L- 
luminosity) for continuum models containing either a thin thermal brems- 
strahlung or a blackbody component and a power-law continuum which is 
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modified by a high energy cutoff factor of the form exp f-(E-E c )/Ep] • 

The variation of a, E c , and Ep across the pulse are now familiar properties 
of other pulsar spectra Including Her X-l (Pravdo et al. 1978) and 
4U0115+63 (Rose et al. 1979). There Is both a spectrally hard pulse region 
and a soft pulse region, as well as a less dramatic spectral evolution 
within each region (Tables 1 and 2). A striking difference, however, between 
4U1 626-67 and Her X-l Is that the spectrally hard section of the pulse has 
a duty cycle of 0.60 In the former compared to only 0.16 in the latter. 

The power emitted in the possible low energy continuum components 
remains fairly constant over the pulse. For the blackbody model the 
temperature and luminosity imply a source radius of x 10 5 d kpc cm. 

This radius is on the order of a neutron star size. A thermal source 
with this size would have an average density of 4 x I0^ 9 d kpc cm’ 3 and 
could be ruled out because of the lack of observed low energy absorption 
or absorption edges. A radius greater than * 10 9 cm is indicated from 
the thermal model. 

Another class of models contains broad emission features rather 
than multiple continuum components (Table 2). A single thin bremsstrahlung 
continuum with a broad (FWHM ^ 15 keV) emission line can also describe 
these spectra. The centroid energy of this feature ranges from ^ 19 
keV to * 7 keV while the power in the line goes from ^ 0.65 to ^ 0.33 of 
the total power. If we assume that the lower energy feature is always 
present, then a two-line model marginally improves the fits for spectra 
1-5, 10. The temperature of the continuum remains fairly constant throughout 
the pulse but changes from -v 10 keV in the single-line model to * 5 keV 
in the two-line model. 
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VI. DISCUSSION 

We have found no evidence for any Doppler variation In the pulse 
phase of 4U1 626-67 that might result from orbital motion with periods 
between 1 and 7 hours. Combining this with the upper limit for periods 
between 7 hours and 300 days obtained by Joss et al. Implies one of the 
following models: 1} the pulsar Is an Isolated neutron star accreting 
matter from an Interstellar cloud; 2) It Is part of a binary system 
whose inclination is low; 3) the orbital period is less than 1 hour 
or greater than 300 days. 

Joss et al. have discussed these possibilities In the light of their 
upper limits between 7 hours and 300 days. They discount the first two 
models as unlikely and propose a binary system with a period of the order 
of a few hours containing a neutron star accreting material from either 
a main sequence dwarf or a low mass white dwarf. For a white dwarf companion 
with < 0.2 M , a sini would be * 0.10 It s, which is at our measured limit. 

O X 

However if the companion were a main sequence M dwarf of mass ^ 0.3 M q 
then our new upper limit indicates that it is of low inclination (i.e. 

% 20°). Thus the nature of this pulsar remains very much an open question. 

Since a neutron star becomes a pulsar because of its superstrong 
magnetic field (e.g. Lamb, Pethlck, and Pines 1973) it is natural to 
assume that the phase-dependent spectral characteristics are caused by 
magnetic effects. However it Is difficult to associate the magnetic field 
energy with a given spectral feature. In the case of 4U1626-67 one might 
conclude that B * 3 x 10 gauss because of the promlnant feature at this 
energy. Or, since another feature may exist near 8 keV the field 
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strength could be half this value with the 'v* 19 keV feature the second 
harmonic. Alternatively If the high energy cutoff Is Interpreted as 
cyclotron absorption (Pravdo et al. 1978) this Implies a field strength 
of * 6 x 10 gauss. This ambiguity results because no viable, compre- 
hensive models that we are aware of now exist. Although there has been a 
considerable theoretical effort begun on this problem. In many cases 
important assumptions of the modelling may not have been justified. This 
particular situation probably requires use of the relatlvlstlcally correct, 
strong-field cross sections (with photon E * — B) In a high temperature 
(T > 20 keV), optically thick plasma (Bussard 1978 and references therein). 

It Is Interesting to note that the broad band (2-60 keV) intensity 
of 4U1 626-67 does not appear to pulse. Over this energy range the narrow 
band pulses destructively interfere. The number of photons remains 
constant, while the energy content of the pulse decreases by 30% with a 
duty cycle of 0.40 (see Table 2). To explain this we must devise a 
phase-dependent mechanism which conserves photon number while sub- 
stantially changing photon energy. Compton scattering does both. However, 
the fact that the transition from hard to soft spectra occurs within 0.1 
phase (Figure 4) indicates that the mechanism can not have performed the 
Isotropization which accompanies normal Compton scatterings. Anisotropic 
Compton scattering in a strong magnetic field (Canuto et al. 1971; Boldt 
et al. 1976) may therefore be an important part of the explanation. 

There may be significance in n comparison of this source with Her X-l. 

In the 2-30 keV Her X-l spectrum a single spectral form can adequately 
describe the spectrum over 84% of the pulse phase with only the normalization 
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left as a free parameter (Pravdo et al. 1978). The normalization changes 
by a factor * 3 Indicating a "background" component of emission which Is 
spectrally but not spatially Isotropic. A Compton scattering shell removed 
from the pulsar surface (Basko and Sunyaev 1976; McCray and Lamb 1976) 
with a suitable geometry could provide this emission ami form the overall 
Intensity envelope (Basko 1977). Another component of emission Is 
viewed directly from the stellar surface In the Her X-l hard X-ray "spectral 
pulse" (Pravdo et al. 1978). We suggest a connection between the presence 
of a background component In Her X-1 and the relatively small duty cycle 
of the spectral pulse. The Compton scattering shell could effectively 
shield the line of sight from most of the surface emission region. The 
7-s pulsar has no background component and no broad-band Intensity pulse. 
Perhaps in this case a Compton shell Is either absent or too small to 
significantly shield any of the surface or to be the source of reflected 
X-rays. This shell is also thought to be the source of the Intense soft 
X-ray (T ^ 10® K) flux from Her X-l (Catura and Acton 1975; Shulman et al. 
1975). Therefore, we would not expect 4U1626-67 to be a strong soft X-ray 
source. The spectral evidence Is not yet complete (l.e. for E < 0.7 keV) 
but does tend to support this final point. 

The spectrum of 4U1626-67 Is so complex that at the present time no 
flnn conclusions about the emitting regions can be made. Detailed 
modelling of the spectral changes across the pulse are required for further 
progress. However, these spectra are only the most extreme example of 
common X-ray pulsar phenomena. Recently features have been detected in the 
X-ray spectra of two pulsars. These have been Interpreted as an emission 
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line In Her X-l (Trumper et el. 1978) and an absorption line In 
4U0115-*63 (Wheaton et al. 1978), originating In the superstrong magnetic 
field {> 10 gauss) at the surface of neutron stars. However in these 
analyses the assumption was made that a d.c. emission component exists 
and Is found In the pnlse Intensity minimum. In both cases cited above 
the Intensity mlnlirao of the pulse was used as a background spectrum 
for the Intensity maximum to obtain a "pulsed spectrum", although Trumper 
et al. obtained similar results of lower statistical significance using 
off-source background. Using the same procedure for the spectra shown In 
Figure 5 would result In dramatic non-physical spectral features. Rose et al. 
(1979) have discussed similar effects In the spectrum of 4U0115+63 in which 
a narrow feature can appear In the net spectrum. We conclude that In 
general, the d.c. component can not be eliminated by this method without 
the danger of subtracting one part of the pulsed emission from another 
and thereby creating spurious features. 

We acknowledge useful discussions with F.K, Lamb, D.Q. Lamb, R. 


Lamb, P.L. Lamb, and J. Cash. 
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CONTINUUM MODELS 
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EMISSION FEATURE MODELS 


One Line Model Two Line Model 
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+ 34 2 

2-60 keV luminosity in units of 10 d kpc erg/s. 



REFERENCES 

Basko, M.M. 1977, Astron Zh. 54, 1050 (transl. In Sov Astron 21, 595). 

Basko, M.M. and Sunyaev, R.A. 1975, Astr. Ag. 42, 311. 

Basko, M.M. and Sunyaev, R.A. 1976, M.N.R.A.S. 175 . 395. 

Boldt, F.A. , Holt, S.S., Rothschild, R.E. and Serlemitsos, P.J. 1976, 

Astr. Aj>. 50, 161. 

Brown, R.L. and Gould, R.J. 1970, Phys . Rev Dl_, 2252. 

Bussard, R.W. 1978, Ph.D. Thesis, U. of Md. 

Canuto, V., Lodenquai, J., and Ruderman, J. 1971, Phys . Rev. 03, 2303. 

Catura, R.C. and Acton, L.W. 1975, Ag. J. ( Letters) 202 , L5. 

Fireman, E.L. 1974, AgL. J. 187, 57. 

Gi«.:oni, R. , Murray, S., Gursky, H., Kellogg, E., Schreier, E., 

Matilsky, T., Kock, D., and Tananbaum, H. 1974, Ag. J. Supple . 27 , 37* 

Ilovaisky, S.A., Motch, C., and Chevalier, C. 1978, Astr . and Ag. 70, LI 9. 

Joss, P.C., Avni, Y., and Rappaport, S. 1978, Ag. 221 , 645. 

Lamb, F.K. , Pethick, C.J., and Pines, D. 1973, Ag. J. 184, 271. 

Markert, T.H., Laird, F.N., Winkler, P.F., Clark, G.W., Hearn, P., Li, 

F.K. , Sprott, G.F. , Bradt, H.V., Lewin, W.H.G., and Schnopper, H.W. 
1978, to be published in Ag. J. ( Supple ). 

Mason, K.O. 1977, M.N.R.A.S . 178 , 81p. 

McClintock, J.E., Canezares, C.R., Bradt, H.V., Doxsey, R.E., and 
Jernigan, J.G. 1977, Nature 270 , 320. 

McCray, R. and Lamb, F.K. 1976, Ag, J. ( Letters ) 204 , LI 75. 

Pravdo, S.H., Proceedings of the Symposium on X-ray Astronomy, XXI COSPAR/ 
IAU, Pergamon Press, Oxford, in press, 1979. 

Pravdo, S.H., Eusscro, R.W., Becker, R.H., Boldt, E.A., Holt, S.S., and 
Serlemitsos, P.J. 1978, Ag. J. 225, 988. 



Rappaport, S., Markert, T. t LI, F.K., Clark, G.W., Jernlgan, J.G., and 
McCUntock, J.E. 1977, A^. J. ( Letters) 217 , L29. 

Rose, L.A., Pravdo, S.H., Kaluzlenski, L.J., Marshall, F.E., Holt, S.S., 

Boldt, E.A., Rothschild, R.E., and Serlemitsos, P.J. 1978, Bull. 

1 1 " i 

of AASjO, 506. 

Rose. L.A., Pravdo, S.H., Kaluzienski, L.J., Marshall, F.E., Holt, S.S., 

Boldt, E.A., Rothschild, R.E., and Serlemitsos, P.J. 1979, submitted 
to Ap. J. 

Rothschild, R., Boldt, E., Holt, S., Serlemitsos, P., Garmire, G., Agrawal , 

P., Riegler, G. , Bowyer, S., and Lampton, M. 1978, NASA TM 79574, 
submitted to Sp. Sci. Instrumentation. 

Schulman, S. , Friedman, H., Fritz, G., Henry, R.C., and Yentis, D.J. 

1975, A£. J. ( Letters) 199 , LI 01. 

Trumper, J., Pietsch, W. , Reppis, C., Voges, W., Staubert, R., and Kendziorra, 
E. 1978, A&. J. ( Letters ) 219 , L105. 

Ulmer, M.P. 1975, Ap, J. 196, 827. 

Wheaton, W. , Primwi, F. , Tsiang, E., Howe, S., Campbell, P., Hoffman, J., 
Lewin, W. , Gruber, D. , Noland, P., Rothschild, R. , Pel ling, R. , and 
Peterson, L. 1978, Bull of AAS 10, 446. 


FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 
Figure 4 


Figure 5 


The master template obtained from these HEAO A-2 obser- 
vations. One sigma uncertainties are shown. The 
amplitude has been normalized relative to the mean count 
rate above the background. The effective area of the HED 
in the 30-60 keV band is still of the order of ^ 50 cm^ and 
the pulse was detected in this band. Although the statistics 
are not good quality the pulse profile is similar both in 
structure, amplitude, and phase to that seen in the 14-30 
keV band. 

The new upper limit of allowed values of a x sini and 

P orb obta ^ ned HEAO A-2 superposed upon earlier SAS 3 

measurements (taken from Joss et al.). The solid and dashed 

lines are loci of the indicated values of the mass and 

radius of the companion for M = 1.5 M and M = 0.5 M 

x ox o 

respectively. Both assume sini = 1. 

The phase-averaged inferred incident spectrum of 4U1 626-67. 
The spectrum of 4U1 626-67 as a function of pulse phase. 

Each spectrum corresponds to a successive 0.1 pulse phase 
starting in the upper left column and continuing down 
first the left and then the right column. The first 
spectrum begins at phase 0 as defined in the light 
curves in Fig. 1. Phase 0 corresponds to JD 2443596.5. 
Superposition of spectra number 3 and 8. 
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